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Innate immune response is the first line of antiviral defense resulting, in most cases, in pathogen clearance with minimal clinical 
consequences. Viruses have developed diverse strategies to subvert host defense mechanisms and increase their survival. In the 
transmissible gastroenteritis virus (TGEV) as a model, we previously reported that accessory gene 7 counteracts the host antivi¬ 
ral response by associating with the catalytic subunit of protein phosphatase 1 (PPlc). In the present work, the effect of the ab¬ 
sence of gene 7 on the host cell, during infection, was further analyzed by transcriptomic analysis. The pattern of gene expression 
of cells infected with a recombinant mutant TGEV, lacking gene 7 expression (rTGEV-A7), was compared to that of cells infected 
with the parental virus (rTGEV-wt). Genes involved in the immune response, the interferon response, and inflammation were 
upregulated during TGEV infection in the absence of gene 7. An exacerbated innate immune response during infection with 
rTGEV-A7 virus was observed both in vitro and in vivo. An increase in macrophage recruitment and activation in lung tissues 
infected with rTGEV-A7 virus was observed compared to cells infected with the parental virus. In summary, the absence of pro¬ 
tein 7 both in vitro and in vivo led to increased proinflammatory responses and acute tissue damage after infection. In a porcine 
animal model, which is immunologically similar to humans, we present a novel example of how viral proteins counteract host 
antiviral pathways to determine the infection outcome and pathogenesis. 


T he order Nidovirales comprises enveloped single-stranded, 
positive-sense RNA viruses and includes the Coronaviridae 
family, which comprises viruses with the largest known RNA ge¬ 
nome (~30 kb) (1, 2). Coronaviruses (CoVs) have been classified 
into three genera— Alphacoronavirus, Betacoronavirus, and Gam- 
macoronavirus (3) —and a fourth, recently proposed, Deltacoro- 
navirus genus (3, 4). These viruses are the causative agents of a 
variety of human and animal diseases. In humans, CoVs produce 
respiratory tract infections, ranging from the common cold to 
severe pneumonia and acute respiratory distress syndrome 
(ARDS) that may result in death (5-9). In animals, CoVs also 
cause life-threatening diseases, such as severe enteric and respira¬ 
tory tract infections, and are economically important pathogens 
(10). However, there is only limited information on the molecular 
mechanisms governing CoV virulence and pathogenesis. 

The 5' two-thirds of the CoV genome encode the replicase 
proteins that are expressed from two overlapping open reading 
frames (ORFs) la and lb (11). The 3' third of the genome contains 
the genes encoding structural proteins and a set of accessory genes, 
whose sequence and number differ between the different species 
of CoV (1, 3). Generally, CoV accessory genes have been related 
with virulence modulation (12). Severe and acute respiratory syn¬ 
drome (SARS)-CoV contains the largest number of accessory 
genes, and it has been proposed that these genes could be respon¬ 
sible for its high virulence (13, 14). A role for some structural 
genes, such as SARS-CoV genes E and 6, on CoV pathogenesis and 
virulence has also been demonstrated (14-18). Nevertheless, in 
general, the function of accessory genes during CoV infection re¬ 
quires further studies (13, 14). 

Double-stranded RNA (dsRNA), produced by RNA viruses as 
a replication intermediate, is a pathogen-associated molecular 
pattern that mediates the activation of well-characterized antiviral 
mechanisms leading to protein synthesis shut down and the stim¬ 
ulation of host innate immunity for initial detection of pathogens 
and subsequent activation of adaptive immunity (19). The path¬ 
way that leads to a block in protein synthesis includes the activa¬ 


tion of double-stranded RNA-dependent protein kinase (PKR), 
leading to eukaryotic translation initiation factor 2 (eIF2a) phos¬ 
phorylation, and the activation of the 2'-5'-oligoadenylate syn¬ 
thetase (2'-5'OAS) and its effector enzyme, the RNase L (RNase 
L), responsible for RNA degradation (19, 20). The host immune 
response triggered by dsRNA is a key component of the innate 
immunity and involves activation of both proinflammatory cyto¬ 
kines and the type I interferon (IFN) system (21, 22). 

There are three main cellular receptors for the detection of 
dsRNA: Toll-like receptor 3 (TLR3), retinoic acid-inducible gene 
I (RIG-I), and melanoma differentiation-associated gene 5 
(MDA5) (22). TLR3 is located in the endosomal membrane of 
antigen-presenting cells, while the cytoplasmic sensors RIG-I and 
MDA5 are the main receptors for viral dsRNA in most cell types 
(20). Recently, degradation of host RNA by RNase L was proposed 
to be an amplifier of the innate immune response by increasing the 
amount of ligand involved in RIG-I and MDA5 recognition (23, 
24). The signaling pathways activated by RIG-I or MDA5 recog¬ 
nition of dsRNA mainly lead to the activation of transcription 
factors IRF3/7 and NF-kB that induce the expression oftype I IFN 
and proinflammatory cytokines (25). This innate immune re¬ 
sponse must be tightly regulated, since there is only a fine line 
separating the induction of a protective antiviral response and an 
exaggerated inflammatory response that can lead to immunopa- 
thology (26). 

Due to the deleterious effects of this response on virus sur- 
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vival, many viruses have developed different strategies that 
counteract the host antiviral responses triggered by dsRNA 
(27). Many of the virus-encoded proteins with this activity 
identified to date interfere with multiple steps of the innate 
response. In addition, some viruses encode more than one gene 
modulating innate immunity (27). CoVs are not an exception 
and encode several proteins affecting type IIFN and proinflam- 
matory cytokines production. Structural proteins, such as nu- 
cleocapsid (N) protein from several CoVs, or SARS-CoV mem¬ 
brane (M) protein have IFN antagonist activity (28-31). The 
modulation of innate immune response by CoV nonstructural 
protein 1 (nspl), nsp3, and nspl6 has also been described. 
Nspl acts by promoting RNA degradation and host proteins 
synthesis suppression (32, 33), reducing both IFN production 
and signaling (34, 35). The antagonist effect of nsp3 is con¬ 
served in different CoV genera and affects IFN and proinflam- 
matory cytokine production, although the mechanism of nsp3 
action has not been determined in all cases (36-39). The IFN 
antagonist effect of nspl6 was recently described, involving a 
mechanism mediated by MDA5 recognition of non-self RNA 
(40). As described above, CoV accessory genes have also been 
related to virulence modulation. Therefore, it could be ex¬ 
pected that some of these genes have a role in innate immunity. 
To date, mouse hepatitis virus (MHV) ns2, 5a, and SARS-CoV 
3b and 6 proteins have been reported as IFN antagonists (28, 
41, 42). Although in general the mechanisms used by accessory 
genes to interfere with the IFN response are not well character¬ 
ized, SARS-CoV protein 6 has been studied in detail. This viral 
protein antagonizes both IFN production (28) and signaling by 
inhibiting signal transduction and activator of transcription 1 
(STAT1) translocation to the nucleus (43). Further, it was re¬ 
cently reported that MHV ns2 protein acts as a 2'-5'-phos¬ 
phodiesterase that reduces the amount of 2'-5'-oligoadeny- 
lates, avoiding the activation of RNase L and, as a consequence, 
reducing RNA degradation during viral infection and type I 
IFN production (24). 

Transmissible gastroenteritis virus (TGEV) is an Alphacorona- 
virus that contains three accessory genes: 3a, 3b, and 7 (44-46). 
TGEV gene 7 is located at the 3' end of the genome and is the last 
ORF. We have recently demonstrated that TGEV protein 7 coun¬ 
teracts host antiviral response and influences virus pathogenesis 
(47). TGEV protein 7 reduces both eIF2a phosphorylation and 
cellular RNA degradation by RNase L (47). The mechanism of 
TGEV protein 7 action is dependent on its binding to cellular 
protein phosphatase 1 (PP1) (47). In addition, infection with a 
mutant virus lacking gene 7 expression (rTGEV-A7) results in 
increased pathological damage compared to the parental (rTGEV- 
wt) virus (47). 

In this work, to understand the molecular mechanisms leading 
to the increased rTGEV-A7 pathogenesis, the role of protein 7 on 
the host cell has been further analyzed by studying differential 
patterns of gene expression during infection with either the wild- 
type or mutant virus. An enhanced proinflammatory response 
was observed in the absence of protein 7, both in vitro and in vivo. 
This increased proinflammatory gene expression was associated 
with elevated levels of macrophage recruitment and activation in 
the infected tissues, being, at least in part, the cause for the en¬ 
hanced tissue damage caused by viral infection in the absence of 
protein 7. 


MATERIALS AND METHODS 

Ethics statement. All animal samples used in the present study were de¬ 
rived from previously published in vivo experiments (47) . As stated in the 
previous publication, these experiments were performed in strict accor¬ 
dance with EU (2010/63/UE) and Spanish (RD 1201/2005 and 32/2007) 
guidelines. 

Virus and cells. Swine testis (ST) cells were grown in Dulbecco mod¬ 
ified Eagle medium supplemented with 10% fetal bovine serum (48). Re¬ 
combinant TGEVs, both parental (rTGEV-wt) and mutant strains lacking 
gene 7 expression (rTGEV-A7), were grown and titrated as previously 
described (47). 

Microarray analysis. One day after achieving confluence, ST cells, 
grown on 35-mm-diameter plates, were mock infected or infected at an 
MOI of 5 with rTGEV-wt or rTGEV-A7. To decrease sample variability, 
nine independent infections were performed in each case. Samples of 
culture supernatants were collected for virus titration as previously de¬ 
scribed (47). Total RNA was extracted, at 6 and 12 hpi, using an RNeasy 
minikit (Qiagen) according to the manufacturer’s instructions. RNAs 
were pooled three by three, obtaining three biological replicates for each 
experimental condition, and RNA integrity was measured in a bioanalyzer 
(Agilent Technologies). RNAs were biotin labeled using a One Cycle Tar¬ 
get labeling kit (Affymetrix). Briefly, cDNA was synthesized from 5 p,g of 
total RNA using an oligo(dT) primer with a T7 RNA polymerase pro¬ 
moter site added to the 5' end. After a second strand synthesis, in vitro 
transcription was performed using T7 RNA polymerase to produce bio- 
tin-labeled cRNA. cRNA preparations (15 p.g) were fragmented at 94°C 
for 35 min into sections of 35 to 200 bases in length and added to a 
hybridization solution (100 mM 4-morpholinopropanosulfonate acid, 1 
M Na + , 20 mM EDTA, 0.01% Tween 20). The cRNAs (10 p,g) were hy¬ 
bridized to Genechip Porcine Genome Arrays (Affymetrix) at 45°C for 16 
h. The arrays were stained with streptavidin-phycoerythrin and read at 
1.56 p.m in a GeneChip Scanner 3000 7G System (Affymetrix). Indepen¬ 
dent microarrays were hybridized for each sample. 

Microarray data analysis. A robust multi-array analysis (RMA) algo¬ 
rithm was used for background correction, normalization, and presenta¬ 
tion of the expression levels (49). Analysis of differential expression was 
performed using a Bayes moderated t test (limma) (50). P values were 
corrected for multiple testing using the Benjamini-Hochberg method 
(false discovery rate [FDR]) (51). The bioconductor packages “affy” and 
“limma” (www.bioconductor.org) were used for these calculations. Due 
to the poor Affymetrix porcine genome array annotation, an alternative 
one based on porcine and human gene homology was used (52). Genes 
were considered differentially expressed when FDR < 0.05. In addition, 
only genes with a fold change of >2 or of <—2 were considered for further 
analysis. 

Functional analysis of microarray results. To understand the biolog¬ 
ical significance underlying the gene expression data, further analysis was 
performed. DAVID (53, 54) was used for the analysis of enriched Gene 
Ontology (55) “biological process” terms in the groups of up- or down- 
regulated genes. 

RNA analysis by quantitative reverse transcription-PCR (RT- 
qPCR). One day after confluence, ST cells, grown on 35-mm-diameter 
plates, were infected at a multiplicity of infection (MOI) of 5. The total 
intracellular RNA was extracted at different hours postinfection (hpi) 
using the RNeasy minikit (Qiagen), according to the manufacturer’s in¬ 
structions. Viral genomic RNA (gRNA) was evaluated by RT-qPCR anal¬ 
ysis, using a custom TaqMan assay and following standard procedures set 
up in our laboratory (56). Cellular gene expression was analyzed using 
TaqMan gene expression assays (Applied Biosystems) specific for porcine 
genes (Table 1). The P-glucuronidase (GUSB) gene was selected as a ref¬ 
erence (housekeeping) gene since its expression remains constant in in¬ 
fected cells (independently of the rTGEV used) compared to noninfected 
cells. Therefore, the expression levels of each gene were corrected by the 
amount of housekeeping gene in each condition. The data were acquired 
with an ABI Prism 7000 sequence detection system and analyzed with ABI 
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TABLE 1 TaqMan assays 


Gene“ 

Assay ID* 1 

Description 

CCL2 (MCP1) 

Ss03394377_ml 

C-C motif chemokine 2 

CCL4 (MIP1B) 

Ss03381395_ul 

C-C motif chemokine 4 

CCL5 (RANTES) 

Ss03648940_ml 

C-C motif chemokine 5 

CXCL9 

Ss03390033_ml 

C-X-C motif chemokine 9 

CXCL11 

Ss03648935_gl 

C-X-C motif chemokine 11 

DDX58 (RIG-I) 

Ss03381552_ul 

DEAD box protein 58 

GUSB 

Ss03387751_ul 

(3 - Glucuronidase 

IFNB1 

Ss03378485_ul 

Interferon beta 

IL-15 

Ss03394854_ml 

Interleukin 15 

IRF1 

Ss03388785_ml 

Interferon regulatory factor 1 

IAK2 

Ss03394066_ml 

Janus kinase 2 

STAT1 

Ss03392296_ml 

Signal transducer and activator of 
transcription 1 

STAT5A 

Ss03394621_ml 

Signal transducer and activator of 
transcription 5A 

TGFB1 

Ss03382325_ul 

Transforming growth factor (31 

TNF 

Ss03391318_gl 

Tumor necrosis factor 

TNFRSF5 (CD40) 

Ss03394339_ml 

TNF receptor superfamily member 5 

VC AMI 

Ss03390912_ml 

Vascular cell adhesion protein 1 


a Alternative gene names are indicated in parentheses. 

b TaqMan assays were used to measure porcine cellular gene expression by RT-qPCR. 
ID, identifier. The “Ss” in these terms refers to Sms scrofa. 


Prism 7000 SDS version 1.2.3 software (Applied Biosystems). All experi¬ 
ments and data analysis were MIQE compliant (57). 

PBMC isolation and cytokine analysis. Lung, blood, and serum sam¬ 
ples were obtained from infected animals (47). Two- to three-day-old 
non-colostrum-deprived piglets, born from TGEV seronegative sows, 
were inoculated with respiratory tropism recombinant viruses (10 7 PFU/ 
pig) according to standard procedures (58). Briefly, animals were infected 
by two different routes (oral and intranasal) in combination. Infected 
animals were monitored daily to detect symptoms of disease and death. At 
0.5, 1, 2, 3, 4, and 5 days postinfection (dpi) two animals per group were 
sacrificed, and the samples were collected. Porcine peripheral blood 
mononuclear cells (PBMCs) were isolated from 8 ml of fresh heparinized 
venous blood by density gradient centrifugation through a Histopaque- 
1077 (Sigma-Aldrich) gradient according to the manufacturer’s recom¬ 
mendations and kept frozen until their use. For their restimulation in 
vitro , PBMCs were cultured (10 6 cells/well) in 24-well plates using RPMI 
medium supplemented with 20% heat-inactivated fetal calf serum (FCS), 
2 mM glutamine, 100 p,M nonessential amino acids, and 100 U of peni- 
cillin-streptomycin/ml. After 18 h, purified TGEV (10 p,g/ml) was added 
for stimulation, and 24 h after stimulation, total RNA was extracted. RNA 
extraction and RT-qPCR measurement of cellular genes was performed as 
described above. 

Cytokine expression analysis from lung samples. At different times 
postinfection, lung sections from infected animals (47) were collected and 
stabilized with RNAlater stabilization reagent (Life Technologies) accord¬ 
ing to the manufacturer’s instructions. Total RNA was extracted by using 
an RNeasy minikit (Qiagen) according to the manufacturer’s instruc¬ 
tions. RT-qPCR measurement of viral gRNA and cellular mRNAs was 
performed as described above. 

Enzyme-linked immunosorbent assay (ELISA). ST cells were in¬ 
fected at an MOI of 5, and the culture medium was harvested at 16 hpi. 
Porcine tumor necrosis factor (TNF), CCL2 and IFN-P protein levels 
were estimated using Swine TNF ELISA kit (Invitrogen), Swine CCL2 
VetSet ELISA development kit (Kingfisher Biotech), and a Pig IFNB 
ELISA kit (Cusabio) according to the manufacturer’s instructions. The 
data were collected from three independent infections. In addition, CCL2 
and TNF were also measured in serum samples according to the manu¬ 
facturer’s recommendations. 


Immunofluorescence of fixed tissues. Representative sections of lung 
tissue were fixed with 4% paraformaldehyde and stored in 70% ethanol at 
4°C (47). Paraffin embedding and sectioning were performed by the his¬ 
tology service in the National Center of Biotechnology (CNB-CSIC; 
Spain). Sections (4 p,m) were deparaffined at 60°C and rehydrated by 
successive incubations in 100% xylol, 100% ethanol, and 96% ethanol. To 
unmask the antigens, tissue sections were boiled in citrate buffer (8.2 mM 
sodium citrate; 1.8 mM citric acid; pH 6.5) for 5 min for macrophage 
detection. Antigen unmasking for granulocyte detection was performed 
by sample incubation with 0.05% trypsin-EDTA for 30 min at 37°C. In all 
cases, after antigen unmasking, the samples were permeabilized with 
0.25% Triton X-100 in phosphate-buffered saline (PBS) for 15 min and 
then blocked with 10% bovine serum albumin (BSA) and 0.25% Triton 
X-100 in PBS for 30 min. The samples were incubated with a monoclonal 
antibody that, in porcine lung samples, is specific for macrophages (4E9 
[kindly provided by J. Dominguez], 1:100) (59), a monoclonal antibody 
that, in porcine samples (60), is specific for granulocytes (Mac387 [Dako], 
undiluted), or a rabbit polyclonal antibody specific for TGEV (1:100). 
Bound primary antibodies were detected with Alexa Fluor 488-conju- 
gated antibody specific for mouse (Invitrogen, 1:250) or Alexa Fluor 594- 
conjugated antibody specific for rabbit (Invitrogen, 1:250), respectively. 
Cell nuclei were stained with DAPI (4',6 , -diamidino-2-phenylindole 
[Sigma], 1:200). Tissues were mounted in Prolong Gold antifade reagent 
(Invitrogen) and analyzed with a confocal fluorescence microscope (TCS 
SP5; Leica). Quantitative analysis of samples was performed with Meta- 
Morph software. Positive pixels for leukocytes were calculated relative to 
positive pixels for TGEV infection. 

Statistic analysis. Two-tailed, unpaired Student t tests were used to 
analyze difference in mean values between groups. All results were ex¬ 
pressed as means ± the standard deviations of the means. P values of 
<0.05 were considered significant. 

RESULTS 

Effect of protein 7 absence on host gene expression. To analyze 
the impact of TGEV protein 7 on host gene expression during infec¬ 
tion, the transcriptomes of rTGEV-wt- and rTGEV-A7-infected cells 
were compared, using porcine microarrays, at two points postinfec¬ 
tion. It is worth noting that under the experimental conditions used 
here (see Materials and Methods), both for the rTGEV-wt and 
rTGEV-A7 viruses, 99% of the cells were infected, and no differences 
in virus titers were observed, as expected (47). MIAME-compliant 
results of the microarrays have been deposited in the Gene Expression 
Omnibus database (GEO [National Center for Biotechnology Infor¬ 
mation] , accession code GSE41756). To select genes expressed at sig¬ 
nificantly different levels in cells infected with either the wild-type or 
mutant viruses, the threshold was established at a fold change of >2 
or <—2 and a false discovery rate (FDR) of <0.05 (Fig. 1A). A com¬ 
parison of infected to noninfected cells, at late times postinfection, 
showed more than 1,500 and 2,500 genes whose expression varied 
significantly between rTGEV-wt and rTGEV-A7 infection (Table 2). 
In contrast, the number of differentially expressed genes was signifi¬ 
cantly reduced when the transcriptomes of cells infected with 
TGEV-wt and TGEV-A7 were compared (Table 2). In relation to 
rTGEV-wt-infected cells, 26 and 301 genes were upregulated, and 84 
and 466 genes were downregulated at 6 and 12 h postinfection (hpi), 
respectively, in rTGEV-A7-infected cells. More than 60% of the dif¬ 
ferentially expressed genes at 6 hpi also changed their expression at 12 
hpi. These results indicated that differential gene activation during 
rTGEV-A7 infection, compared to that induced by rTGEV-wt, was 
maintained at 12 hpi, although additional differentially expressed 
genes were noted at later times postinfection. 

Genes differentially expressed during TGEV infection (both 
up- and downregulated), in the presence or absence of protein 7, 
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DNA replication 
Phosphorylation 
Transition metal ion transport 
Regulation of hydrolase activity 
RNA modification 

Negative regulation of macromolecule metabolic process 
Carbohydrate biosynthetic process 
Organelle fission 
Response to organic substance 
Regulation of cell killing 
Nuclear division 
Mitosis 

Regulation of transcription from RNA polymerase II promoter 
Positive regulation of nuclear division 
Chronic inflammatory response 
Cellular macromolecular complex subunit organization 
Regulation of myeloid cell differentiation 
Negative regulation of gene expression 
Positive regulation of nitrogen compound metabolic process 

Cell proliferation 
Positive regulation of developmental process 
Nucleosome assembly 
Regulation of innate immune response 
Leukocyte adhesion 
Posttranscriptional regulation of gene expression 
Cytokine-mediated signaling pathway 
Leukocyte activation 
Regulation of programmed cell death 
Response to bacterium 
Unfolded protein response 
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Inflammatory response 
Ribosome biogenesis 
Immune response 
Defense response 
Response to virus 

io - 2 10- 3 10- 4 10- 5 10* 10- 7 10* 10* 



P-VALUE 

FIG 1 Effect of TGEV protein 7 on host gene expression. (A) Comparison of gene expression in ST cells at 6 and 12 hpi using microarrays. The data visualization and 
threshold set up for results filtering were performed with FIESTA viewer (J. C. Oliveros, 2007, http://bioinfogp.cnb.csic.es/tools/FIESTA/index.php). The graphs 
represent the normalized ratios for rTGEV-wt-infected versus mock-infected cells, rTGEV-A7-infected versus mock-infected cells, and rTGEV-A7-infected versus 
rTGEV-wt-infected cells. Dark gray spots indicate upregulated probes (fold change, >2), while medium gray spots indicate downregulated ones (fold change, <—2). 
Only genes with an FDR of <0.05 were considered as candidate genes. (B) Candidate genes that were differentially expressed in rTGEV-A7-infected cells compared to 
rTGEV-wt-infected ones, at 12 hpi, were grouped with reference to their GO biological process terms. Numbers on the x axis indicate DAVID FDR values. 
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TABLE 2 Analysis of host gene expression using porcine microarrays 


Comparison 

Time point 
(hpi) 

No. of probes* 


Upregulated 

Downregulated 

WT vs Mock 

6 

373 

148 


12 

1,119 

768 

A7 vs Mock 

6 

289 

127 


12 

1,259 

1,782 

A7 vs WT 

6 

31 (26) 

98 (84) 


12 

346 (301) 

525 (466) 


a Numbers in parentheses represent the numbers of different genes, since one gene 
could be present in the array with more than one probe. Upregulated, fold change > 2, 
FDR < 0.05; downregulated, fold change <—2, FDR < 0.05. 


were grouped, using DAVID software (53, 54), according to the 
biological processes in which they could be potentially involved. 
Most of the genes were associated with responses against viruses, 
host defense response, and immune response (Fig. IB) and there 
was significant overlap among the genes of these three groups. The 
genes included in response to virus gene ontology (GO) group 
mainly affected the immune response, the IFN response, and in¬ 
flammation, and most of them were significantly upregulated dur¬ 
ing rTGEV-A7 infection compared to that by rTGEV-wt (Fig. 2A). 
Interestingly, analysis of the expression patterns during rTGEV-wt 
and rTGEV-A7 infections showed that, in general, genes activated 
during rTGEV-wt infection were also activated during rTGEV-A7 
infection, although earlier or to a higher level than in rTGEV-wt in- 


IMMUNERESPONSE 

+ 6 / 7 ^ 
+4.5 
+4.5 
+4.3 
+3.9 
+2.5 
+ 2.2 
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VESICLE TRAFFIC 
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STAT5A +2.3 
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IFN RESPONSE 


ZlFNp 

+24.^\ 

DDX58 (RIG-1) 

+5.4 

IFIH1 (MDA5) 

+3.2 

ISG20 

+3.2 

G1P3 

+2.9 

IRF1 

+2.4 

STAT1 

+2.2 

MX2 

+2.0 

G1P2 

+2.0 

\IRF7 

+2.oy 

COMPLEMENT ACTIVATION 

BF 

+2.2 


NUCLEAR TRANSPORT 

KPNA2 -2.3 

RANBP5 -2.3 


INFLAMMATION 


/CCL2 (MCP-1) 

+7^9\ 

TNF 

+6.8 

CCL4 (MIP-lp) 

+5.1 

MAPK9 (JNK2) 

+4.1 

CCL19 

+3.5 

CXCL9 

+3.2 

VCAM1 

+2.7 

ICAM1 

+2.6 

CCL5 (RANTES) 

+2.6 

NFkBIA 

+2.5 

OLR1 

+2.4 

IL8 

+2.0 

MCP 

-2.0 

1 ILIA 

-2.y 

APOPTOSIS 

c 


TNFRSF11B 

+2.2 

TNFRSF21 

-2.2 

TMEM23 

-3.5 


V_/ 


CXCL11 

IL15 

PTPRC 

HLA-A 

TNFRSF5 (CD40) 

SOCS1 

ZC3HAV 

ICOSL 

TAPI 

LYZ1 

IL1RN 

CXCL6 



FIG 2 Host genes differentially expressed in rTGEV-A7 versus rTGEV-wt-infected cells. (A) Candidate genes included in response to virus GO group were 
classified according to their main biological functions. Black and gray lettering is used to indicate up- and downregulated genes, respectively. The numbers 
indicate the fold change for each gene. For genes recognized with more than one probe, the value corresponding to the greatest up- or downregulation is 
represented. (B) Thirteen candidate gene mRNAs were evaluated by RT-qPCR using specific TaqMan assays. RNA was extracted from rTGEV-wt- or rTGEV- 
A7-infected cells at 12 hpi. Gray bars represent the fold change obtained from microarray data. Black bars represent the data obtained by RT-qPCR. GUSB mRNA 
levels were used as an endogenous control in all cases. Error bars indicate the standard deviations from three independent experiments. Imm. Resp., immune 
response; Trans., transcription. 
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TABLE 3 Candidate genes selected for further study* 


Immune category 

Cytokine 

Function 

Fold change 

Innate immunity 

IFN-0 

Antiviral, increased MHC class I expression 

+ 24.3 

Thl response 

IL-15 

Stimulates growth and proliferation of T and NK cells 

+ 4.5 


TNF 

Local inflammation, endothelial activation; stimulates cell proliferation and induces cell differentiation 

+ 6.8 


CXCL9 

Chemotactic for activated T cells; affects the growth, movement, or activation state of cells 
participating in immune and inflammatory responses 

+ 3.2 


CXCL11 

Chemotactic for interleukin-activated T cells; role in diseases involving T-cell recruitment 

+ 6.7 

Chemokine 

CCL2 

Attracts monocytes and basophils; related to diseases characterized by monocytic infiltrates 

+ 7.9 


CCL4 

Inflammatory and chemokinetic properties 

+5.1 


CCL5 

Attracts monocytes, memory T-helper cells, and eosinophils; causes the release of histamine from 
basophils and activates eosinophils 

+2.6 

Others 

VCAM1 

Important in cell-cell recognition; influences leukocyte-endothelial cell adhesion 

+ 2.7 

Negative controls 

TGFB1 

Controls proliferation, differentiation, and other functions in many cell types 

-1.15 


JAK2 

Mediates essential signaling events in both innate and adaptive immunity 

+ 1.67 


a The selection was based on their differential expression levels and porcine TaqMan assay availability. 


fection (data not shown). Several genes from each group were se¬ 
lected for RT-qPCR analysis, taking into account the availability of 
TaqMan assays for porcine genes detection. In all cases, RT-qPCR 
results confirmed the microarray results, validating the effects ob¬ 
served in the microarray analysis (Fig. 2B). Altogether, these data 
suggested that host innate immune response triggered by TGEV in¬ 
fection was magnified in the absence of protein 7. 

Effect of protein 7 absence on innate immunity gene expres¬ 
sion. We have previously shown that, due to protein 7 absence, 
rTGEV-A7 infection promoted an intensified host dsRNA-in- 
duced antiviral response (47). Taking this observation into ac¬ 
count, a set of nine differentially expressed candidate genes in¬ 
volved in innate immunity (Table 3) was selected to study their 
expression kinetics by RT-qPCR. To discard a general transcrip¬ 
tion upregulation during rTGEV-A7 infection, which could affect 
the expression of all cellular genes, the mRNA levels of two other 
genes involved in immune response, such as JAK2 and transform¬ 
ing growth factor (3 (TGF-(3), were also evaluated as controls. No 
differences in virus titers (data not shown) or viral gRNA accumu¬ 
lation (Fig. 3) were observed between the parental and mutant 
viruses, as expected (47). All of the selected candidate genes in¬ 
creased their expression during TGEV infection (Fig. 3). Interest¬ 
ingly, all of them were significantly upregulated in cells infected 
with rTGEV-A7, in comparison to cells infected with rTGEV-wt 
(Fig. 3). In contrast, control genes JAK2 and TGF-p were not 
differentially expressed in the absence or presence of TGEV pro¬ 
tein 7 (Fig. 3), indicating that the increased expression of candi¬ 
date genes during rTGEV-A7 infection was not the consequence 
of a general induction of host transcription. 

To confirm that the elevated mRNA levels correlated with an 
increased protein accumulation, and attending to the availability 
of assays for the specific detection of porcine proteins, the levels of 
three cytokines (TNF, CCL2, and IFN-(3) were measured during 
TGEV infection. TNF, CCL2, and IFN-p levels were quantified in 
the medium of infected ST cells using commercial ELISAs. Cells 
infected with rTGEV-A7 virus secreted higher levels (Fig. 4) of 
TNF (upper panels), CCL2 (medium panels) and IFN-(3 (lower 
panels) compared to rTGEV-wt-infected cells. Altogether, the 
data confirmed the correlation between mRNA and protein ex¬ 
pression levels of innate immunity genes induced by rTGEV in¬ 
fection. 


In vivo expression of proinflammatory cytokines in the ab¬ 
sence of TGEV protein 7. It was previously shown that lung dam¬ 
age in rTGEV-A7-infected pigs was greater than that observed in 
rTGEV-wt-infected animals (47). The observed lesions could be 
correlated with acute inflammation in the infected tissue (26). As 
described above, the data obtained in tissue culture indicate that 
rTGEV-A7 infection led to an increased expression of innate im¬ 
munity-related genes, in particular of proinflammatory cytokines. 
To analyze whether this was also the case in vivo, cytokine expres¬ 
sion in PBMCs from infected animals was studied by RT-qPCR. 
After in vitro restimulation of PBMCs from rTGEV-wt- and 
rTGEV-A7-infected animals with purified TGEV, an increase in 
the expression of several proinflammatory cytokines, such as 
TNF, CCL2, CCL5, and CXCL11, was observed (Fig. 5). TGEV 
infection of PBMCs was not observed, in contrast to what has been 
previously described, indicating that TGEV could infect pDCs 
(61-63). This cell type represents ca. 0.5% of the cells in a PBMC 
preparation (64), and pDC enrichment methods are required for 
the significant detection of TNF and type IIFN production by this 
cell type (62). Therefore, the observed increase in proinflamma¬ 
tory cytokines production was most likely due to re-exposure to 
TGEV. The levels of cytokine mRNAs were higher in PBMCs from 
rTGEV-A7-infected animals than in those from rTGEV-wt-in- 
fected animals (Fig. 5), in agreement with the results obtained in 
tissue culture. These data strongly suggested that rTGEV-A7 virus 
also induced a higher expression of proinflammatory cytokines in 
vivo. 

To further analyze the proinflammatory cytokine levels in the 
infected animals, mRNA levels in lung samples were measured by 
RT-qPCR. Since the infection extent in the different lung sections 
was different, viral gRNA was also quantified, and cytokine mRNA 
levels were related to viral gRNA levels (Fig. 6). At 1 dpi, the 
mRNA levels of TNF, CCL2, and CCL5 were higher in rTGEV- 
A7-infected lungs than in rTGEV-wt-infected lungs (Fig. 6). In 
the case of CCL2, this increase was also found at 2 dpi (Fig. 6, 
middle panel). Moreover, CCL2 and TNF protein levels were 
measured in the sera from infected animals. Sera from rTGEV- 
A7-infected animals contained higher levels of these cytokines 
than those from rTGEV-wt-infected animals, particularly at early 
dpi (Fig. 7). Altogether, the data indicated that rTGEV-A7 virus 
also induced an exacerbated proinflammatory response in vivo. 
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FIG 3 Expression kinetics of host genes involved in inflammation. Nine candidate genes involved in inflammation and innate immunity were selected for 
analysis of their expression kinetics by RT-qPCR. ST cells were infected with rTGEV-wt or rTGEV-A7, and the total RNA was extracted at the indicated time 
points. As a control for rTGEV replication, the viral gRNA levels were also evaluated. TGF-p and JAK2 were selected as negative control genes that did not change 
their expression levels in rTGEV-A7-infected cells compared to rTGEV-wt-infected cells. GUSB mRNA levels were used as an endogenous control in all cases. 
Error bars indicate the standard deviations from three independent experiments. **, P < 0.01; *, P < 0.05. 
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FIG 4 Host protein and mRNA levels in rTGEV-A7 infection. Quantification of 
porcine TNF (upper panels), CCL2 (middle panels), and IFN-(3 (lower panels) 
accumulation in mock-infected cells (white) or during rTGEV-wt (black) or 
rTGEV-A7 (gray) infections. At 16 hpi, mRNA accumulation was measured by 
RT-qPCR (left), and protein accumulation was measured by ELISA (right). Error 
bars indicate the standard deviations from three independent experiments. ***, 
P < 0.001; **,P< 0.01. 


The expression of proinflammatory cytokines during infection 
could lead to an enhancement of tissue damage due to the activa¬ 
tion and recruitment of leukocytes, particularly macrophages and, 
in some infections, neutrophils (65). These activated cells may 
amplify the damage produced by the pathogen, since they secrete 
proinflammatory cytokines potentially involved in host-mediated 
immunopathology (65,66). To test this possibility, the presence of 
granulocytes and macrophages in the lungs of rTGEV-wt- and 
rTGEV-A7-infected animals at 1 and 4 dpi was evaluated by im¬ 
munofluorescence. At 1 dpi, an increase in granulocytes was ob¬ 
served in both rTGEV-wt- and rTGEV-A7-infected tissues (Fig. 
8A) that correlated with the damage observed by histopathology 
(47). The number of granulocytes seemed to decrease at 4 dpi in 
both infections (Fig. 8A), although the lung damage in rTGEV-A7 
infection was increased compared to the one observed at 1 dpi 
(47). These data suggested that granulocytes were rapidly re¬ 
cruited at the sites of infection, as expected. The contribution of 
macrophages to tissue damage was analyzed by using the antibody 
4E9/11 that specifically labeled porcine lysosomal associated 
membrane protein 1 (LAMP-1, CD107a), which mainly detects 
cells in an active phagocytosis state (59, 67). Therefore, this anti¬ 
body allowed the simultaneous detection of both activated resi¬ 
dent macrophages and circulating monocytes recruited into the 
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FIG 5 Expression of cytokines by PBMCs from infected animals. PBMCs from 
rTGEV-wt (black)- or rTGEV-A7 (gray)-infected animals, extracted at 0.5 and 
3 dpi, were cultured and stimulated with TGEV. At 24 h after stimulation, total 
RNA was extracted, and the cytokine mRNA levels were analyzed by RT-qPCR. 
GUSB mRNA levels were used as endogenous control in all cases. Error bars 
indicate the standard deviations from two different animals. ***, P < 0.001; **, 
P < 0.01; *,P< 0.05. 
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FIG 6 Expression of cytokines in lungs from infected animals. Lung fragments 
from rTGEV-wt (black)- or rTGEV-A7 (gray)-infected animals were collected 
at 1, 2, and 4 dpi. Total RNA was extracted, and the viral gRNA accumulation 
and cytokine mRNA levels were analyzed by RT-qPCR. TNF (upper panel), 
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FIG 7 Cytokine levels in sera from infected animals. Quantification of porcine 
CCL2 (upper panel) and TNF (lower panel) accumulation in sera from mock- 
infected piglets (white) or animals infected with rTGEV-wt (black) or rT- 
GEV-A7 (gray). Cytokine accumulation was measured by ELISA at the indi¬ 
cated times postinfection. Error bars indicate the standard deviations from 
three independent experiments. **, P < 0.01; *,P < 0.05. 


site of infection, providing that they were subsequently activated. 
Other porcine macrophage markers, such as CD 163, CD68 and 
CD 172a, were tested. Unfortunately, the antibodies specific for 
these markers did not work in paraffin-embedded lung sections. 
An increased number of activated macrophages was observed in 
rTGEV-A7-infected lungs at 4 dpi compared to those infected 
with rTGEV-wt (Fig. 8B). Although it was previously reported 
that TGEV can replicate in alveolar macrophages (61), double- 
labeled cells (4E9/11 1 TGEV + ) were not detected in the lung sec¬ 
tions at any time point analyzed. It was previously observed that 
rTGEV-A7 titers in lung were higher than those of the rTGEV-wt 
virus at early times postinfection (47). Therefore, the contribution 
of differences in tissue infection levels to leukocyte recruitment 
and activation was a realistic possibility. To clarify this issue, gran¬ 
ulocyte and macrophage recruitment and activation, during 
rTGEV-wt and rTGEV-A7 infections, were estimated relative to 


CCL2 (middle panel), and CCL5(lower panel) mRNA levels were determined 
relative to the viral gRNA in all cases. In addition, the GUSB mRNA levels were 
used as an endogenous control in all cases. Error bars indicate the standard 
deviations from two different animals. ***, P < 0.001; **, P < 0.01; *, P < 0.05. 
r.u., relative units. 
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FIG 8 Leukocyte detection in lung sections from infected animals. Lung sections from noninfected (Mock) and rTGEV-wt (WT)- or rTGEV-A7 (A7)-infected 
piglets, at 1 and 4 dpi, were labeled with a monoclonal antibody specific for granulocytes (A) or macrophages (B) (red) and a polyclonal antibody specific for 
TGEV (green). Cell nuclei were labeled with DAPI (blue). Pictures were obtained by using confocal microscopy with a X40 objective lens. 



tissue infection levels in each case. Granulocyte recruitment was 
similar in both wild-type and mutant virus-infected tissues 
(Fig. 9). In contrast, at 4 dpi rTGEV-A7-infected tissues showed a 
significant increase in macrophage recruitment and activation 
compared to rTGEV-wt-infected tissues (Fig. 9). Altogether, these 
data suggested that the enhancement of lung damage produced by 
rTGEV-A7 infection was, at least in part, due to the preferential 
recruitment and activation of macrophages, most likely due to the 
increased levels of proinflammatory cytokines such as TNF or 
CCL2. 

DISCUSSION 

We have previously shown that an rTGEV lacking protein 7 ex¬ 
pression was more virulent and caused increased pathology than 
the parental virus (47). To analyze the potential mechanism un¬ 
derlying this enhanced pathology, the patterns of gene expression 
after infection by each of these viruses were analyzed using mi¬ 
croarrays covering the complete porcine genome. A marked up- 
regulation of proinflammatory cytokine mRNA expression was 


observed in infections with the virus deficient in protein 7. The 
identification of elevated TNF, CCL2, and IFN-(3 confirmed the 
increased proinflammatory pattern at the protein level. Further¬ 
more, similar results were obtained in in vivo infections, indicat¬ 
ing that the presence of protein 7 reduced inflammatory changes 
after TGEV infection both in cell cultures and in vivo. 

Viruses are good tools for understanding the molecular mech¬ 
anisms modulating inflammation, especially signaling pathways 
that increase disease severity. The increased inflammation ob¬ 
served after rTGEV-A7 infection, caused by an exacerbated innate 
immune response and leading to an enhanced pathogenesis, was 
also described for human viruses infecting the respiratory tract, 
such as respiratory syncytial virus (RSV) or influenza vims. TGEV 
is a virus with enteric and respiratory tropism. Lung and gut in¬ 
fection by virulent TGEV caused significant inflammation in both 
tissues, and animal death is mainly due to the severe unbalance of 
Na + and K + ions caused by the clinical manifestation of the in¬ 
fection (68). It is important to note that the work described here 
was performed with the cell culture-adapted TGEV used in the 
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FIG 9 Quantification of the leukocyte/infection ratio in lung sections from 
infected animals. The leukocyte signal (immunofluorescence red channel) was 
estimated in lung sections from infected animals. The positive area for granu¬ 
locytes (upper panel) or macrophages (lower panel) signal relative to the pos¬ 
itive area for TGEV infection (immunofluorescence green channel) signal was 
calculated using MetaMorph software. Error bars indicate the standard devi¬ 
ations from 20 observed fields in 10 independent samples. ***, P < 0.001. r.u., 
relative units. 


previous study (47), which only displays respiratory tropism caus¬ 
ing lung damage and no apparent gut infection. TGEV is a porcine 
virus, and the work presented here has been performed using the 
natural host, which is immunologically more similar to humans 
(>80%) than mice (<10%) (69). Therefore, our findings might 
be informative for the conserved pathways leading to increased 
pathology both in pigs and humans. 

TGEV protein 7 reduced the dsRNA triggered antiviral re¬ 
sponse (47). As a consequence, rTGEV-A7 infection caused en¬ 
hanced cell death, cellular RNA degradation, and protein synthe¬ 
sis shutoff (47). The results presented here were most likely not 
conditioned by these effects, since both microarray and RT-qPCR 
analyses were performed using the same amount of total intracel¬ 
lular RNA of equal quality, extracted from living cells attached to 
the plate. In addition, microarray data normalization corrects any 
differences in mRNA amount between samples, and in RT-qPCR 
the mRNA levels always referred to the amount of GUSB that 
correlated with the percentage of living cells. In the case of cyto¬ 
kine measurements in the supernatants of infected cells, cell death 
most likely did not significantly affect quantifications, since simi¬ 
lar differences between rTGEV-wt and rTGEV-A7 viruses were 


obtained when intracellular protein extracts were used for ELISAs 
(data not shown). In addition, up to a 20-fold increase in mRNA 
accumulation led to an increase no more than 3.5-fold in protein 
levels when rTGEV-wt and rTGEV-A7 infections were compared 
(Fig. 4). Similarly, differences up to 10 3 -fold in mRNA levels led to 
no more than 10-fold differences in protein accumulation, when 
comparing mock-infected cells to rTGEV-wt-infected cells (Fig. 
4). This result suggested that shutoff caused by the rTGEV-A7 
virus was not affecting cytokine measurements. 

The infection of swine by rTGEV-A7 virus caused higher lung 
damage than that caused by rTGEV-wt (47). Neutrophils and 
macrophages have been proposed to cause lung damage during 
most cases of acute lung injury and ARDS (70) . Granulocytes were 
recruited to the TGEV-infected lungs, both in the presence and in 
absence of protein 7, suggesting that they play a role in TGEV- 
induced inflammation. In vivo, rTGEV-A7 infection caused an 
increased expression of proinflammatory cytokines such as TNF, 
CCL2, and CCL5, whose main function is macrophage recruit¬ 
ment and activation (26, 71, 72). In addition, at least in vitro, 
rTGEV-A7 infection induced the expression of CD40 receptor 
(Fig. 2B), which has been involved in macrophage activation (73). 
In agreement with these results, in the lungs from rTGEV-A7- 
infected animals, an increase in macrophage recruitment and ac¬ 
tivation was observed compared to the rTGEV-wt-infected ani¬ 
mals. In line with these data, MHV infection of CCL5 receptor 
knockout mice causes a reduced pathology due to a decrease in 
macrophage recruitment (74). In addition, it was recently re¬ 
ported that coronavirus infection of transgenic mice expressing 
CCL2 led to an enhanced pathology leading to death, caused by a 
dysregulated immune response without effective virus clearance 
(75). Similar observations were obtained after influenza virus in¬ 
fection of CCL2 (CCR2~ /_ ) or CCL5 (CCR5~ /_ ) receptors 
knockout mice (76). Lung pathology was reduced in influenza 
virus-infected CCR2~ / ~ mice, since monocyte recruitment is se¬ 
verely impaired in these animals (76). In contrast, CCR5~ /_ in¬ 
fected mice showed increased mortality associated with elevated 
macrophage infiltration in the lungs due to an increased expres¬ 
sion of CCL2 (76). The data obtained in the present work sug¬ 
gested that macrophages were involved in the enhanced inflam¬ 
mation produced during infection in the absence of protein 7 and 
were in agreement with the role of CCL2 in the immunopathology 
mediated by macrophage recruitment. 

An exacerbated proinflammatory cytokine production and an 
excessive immune cell recruitment leading to tissue destruction 
contributing to virus caused pathology, similar to that observed 
after rTGEV-A7 infection, has been described for several viral in¬ 
fections (77-79). This immunopathology, known as “cytokine 
storm,” could be the cause for the extreme virulence of several 
viruses, such as pandemic influenza virus H5N1 or SARS-CoV 
(80). Once started, increased proinflammatory cytokines could 
continue driving immunopathology progression, even in the ab¬ 
sence of viral replication. In fact, lung damage in SARS-CoV- 
infected patients persists after virus titer reduction, suggesting 
that pathology is mainly caused by the immune response (81). 
Similarly, in RSV infections, the severity of the infection has been 
correlated with CCL2 and CCL5 expression (72). In addition, cel¬ 
lular recruitment, mediated by cytokine expression, has been con¬ 
sidered responsible for damaging both infected and uninfected 
areas of the lung (26). In line with these observations, rTGEV-A7 
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infection caused a “cytokine storm” that led to a progression in 
lung damage (47). 

In the absence of protein 7, an increased expression in IFN-p, 
IFN-stimulated genes (ISGs), and proinflammatory genes was ob¬ 
served. Similar upregulation of genes was reported for viruses with 
mutations in virus-encoded IFN antagonists, such as influenza 
virus with the mutated NS1 gene (82, 83). To date, virus-encoded 
IFN antagonists were mainly analyzed in overexpression studies, 
and the activity of only around half of these antagonists was dem¬ 
onstrated in the context ofviral infection (27). The IFN antagonist 
activity of TGEV protein 7 decreasing IFN-p production was 
demonstrated here, in the context of viral infection, although 
more studies are required to further determine the characteristics 
and mechanism of IFN antagonism by protein 7. 

As expected for a virus defective in an IFN antagonism path¬ 
way, the rTGEV-A7 mutant virus was more sensitive to IFN-p- 
induced antiviral effects (data not shown). Nevertheless, the rT- 
GEV-A7 virus underwent efficient replication despite increased 
IFN-P production. This was explained by the amount of IFN-p 
produced by ST cells, being 10 3 -fold lower than the minimal con¬ 
centration required for decreasing TGEV replication (data not 
shown). Most likely, this result was a consequence of the presence 
of other virus-encoded IFN antagonists in both rTGEV-wt and 
rTGEV-A7 genomes. 

We have previously shown that TGEV protein 7 limited RNase 
L activation and eIF2a phosphorylation through its binding to 
PP1 phosphatase (47). RNase L has been involved in the IFN re¬ 
sponse, since cells deficient in this enzyme, infected with different 
viruses, produce lower IFN-p amounts than normal cells (23). It 
was proposed that the RNA degradation products generated by 
RNase L are recognized by RIG-I, acting as amplifiers of IFN pro¬ 
duction (23). In fact, the direct implication of RNase L on IFN 
production has been recently demonstrated using an MHV acces¬ 
sory protein ns2 mutant virus (24) . In addition, PKR has also been 
recently involved in the amplification of innate immune response 
through a translational control mechanism, dependent on eIF2a 
phosphorylation, leading to increased IFN-p production and 
NF-kB activation (84). Therefore, the increased IFN-p and pro- 
inflammatory cytokines production observed during rTGEV-A7 
infection, in the absence of protein 7, could be explained in the 
context of the previously proposed TGEV protein 7 mechanism of 
action, highlighting the role of RNase L and eIF2a phosphoryla¬ 
tion in innate immunity. The results presented here confirm the 
role of CoV accessory genes in the modulation of innate immune 
responses during infection. The tight regulation of deleterious in¬ 
flammatory responses by virus-encoded proteins seems to modu¬ 
late both virus and host survival. 
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